Conductance of a photochromic molecular switch with graphene leads 
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We report a full self-consistent ab initio calculation of the conductance of a diarylethene-based 
molecular switch with two graphene electrodes. Our result show the contributions of the resonant 
states of the molecule, of the electrode density of states, and of graphene unique features such as 
edge states. The conductivities are found to be significantly different for the two photochromic 
isomers at zero and finite applied bias. Further we point out the possibility of causing the switching 
by the application of a large potential difference between the two electrodes. 



Organic electronics has rapidly grown into a funda- 
mental field whose potential is still to be fully developed 
and exploited. Research focuses on novel functional or- 
ganic materials and existing applications already com- 
prise, among others, nanoscale electronic devices such as 
thin film transistors and diodes, solar cells, integrated cir- 
cuits on flexible substrates, carbon nanotube field effect 
transistors and molecular switches. A promising class of 
organic molecular-scale photoswitching units is provided 
by photochromic molecules which are able to switch be- 
tween two chemical structures when irradiated by light 
of appropriate wavelengths. Such states correspond to 
photochemically interconvertible isomers [l[. 

A recently developed class of molecules, named di- 
ary lethencs, displays a photochromic, thermally irre- 
versible, reactivity 0, Most of them exhibit light 
induced switching both in solution and in single crys- 
tals @, 0, El by changing reversibly in aromaticity dur- 
ing the closed-open (open-closed) configuration transi- 
tion, determined by visible and ultraviolet light, respec- 
tively |f|. This property raised much interest and stim- 
ulated research on diarylethene derivatives as suitable 
photochromic molecular switches @, H|. Investigations 
of diarylenthenc junctions with gold electrodes showed 
that conductivity increases by two orders of magnitude 
in the closed isomer compared to that of the open one Q . 
Nonetheless several diarylethene derivatives on gold leads 
only display photochromic switching from the closed to 
the open state after irradiation with visible light but 
not the reverse process by UV [Io|, II owing to the 



spacer reversible light-induced switching resulted both on 
Au(lll) [l3| and gold nanoparticles 14\. Photochromic 
diarylethene derivatives are also investigated on single 
walled carbon nanotube (SWCNT) as perspective single 
molecule devices with a more accurately defined contacts 
to the electrodes [lH . For such system it was found that 
switching from insulating to conducting, i.e. from open 
to closed molecular configuration occurs, but not the re- 
verse, unlike in the case of Au [l6|. The timely charac- 
ter of this topic is confirmed by the wealth of theoreti- 
cal studies on the conductance of photochromic switches 
both with gold HHEjEII and SWCNT electrodes [13]. 



fast quenching of the photo-excited hole of the HOMO 
state of the open isomer into gold |12|. By connect- 



ing the metal anchoring sulfur atom by a phenyl group 



In this paper we investigate the conductance of a di- 
aylethcnc switch sandwiched between two semiinfinitc 
graphene electrodes by the Landauer formalism recast in 
terms of the non-equilibrium Green's function (NEGF) 
within density functional theory (DFT) [23|. In fact, the 
two-dimensional (2D) form of carbon has just revealed 
several intriguing fundamental properties and huge po- 
tential in nanoelectronics. In particular, high chemical 
stability, low resistivity and mechanical strength suggest 
graphene, a more cost-effective material than carbon nan- 
otubes, as an alternative component for electrodes in 
electronic devices [2~ij . Furthermore in graphene elec- 
trons can travel ballisticalfy up to a fim. (25| , an d electron- 
phonon scattering intensity is very weak [26J. This im- 
plies a carrier mobility larger than that of the inorganic 
semiconductor with the largest mobility, i.e. InSb, and 
that of carbon nanotubes. If defects and impurities are 
eliminated, such mobility is about 2 • 10 5 cm 2 /Vs [25| . 
but in practice graphene sheets have irregular shapes 
and contain impurities and defects. To overcome this 
difficulty we point out a promising bottom up experi- 
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mental technique capable to synthetize graphene start- 
ing from small corencne blocks and to build well defined 
molecule electrode contacts [27]]. Transparent, conduc- 
tive, and ultrathin graphene films, have been demon- 
strated to be alternative metal-oxide electrodes for solid- 
state dye-sensitized solar cells [28| . 

Our calculations show that at zero bias the peaks in 
conductance of the diarylethene moiety attached to two 
graphene electrodes reflect not only the molecule elec- 
tronic states broadened into resonances by the coupling 
to the leads, but also the features of graphene edge states. 
Near the Fermi energy (Ep) the conductance of the closed 
isomer is much lower in the energy interval between the 
HOMO and the LUMO resonances of the molecule, fol- 
lowing the linear dependence of the electrode density of 
states (DOS) proportional to \E—Ep\. In the same range 
the conductance of the open isomer is about zero. If we 
apply a finite bias we observe that conductance is allowed 
within different energy windows for the open and closed 
isomers. 

We studied the open and closed structures of two dif- 
ferent diarylethene molecules, shown in Fig. [TJ The cen- 
tral dithicnylcthene part is the same for each molecule, 
and it is attached to the electrodes through — R— with 
R=CH 2 (we call this Al) and R=CO-NH (referred to 
as A2). We also considered longer isomers, which have 
two phenyl linkers and are connected in the same fashion, 
namely Bf and B2. We started with optimized structures 
of the molecules with the two R groups connected to H 
atoms. Then the two H atoms were removed and the 
molecules were embedded between two graphene sheets. 
The graphene electrodes were taken to have zig-zag edges 
with the valencies of the edge C atoms satisfied by H 
atoms. We did this because C atoms at zig-zag termina- 
tions display stronger affinity towards radicals, and hence 
form stronger bonds with H atoms. Graphene with zig- 
zag edges have unique edge state near the Fermi level 
[29j . Further, we also optimized the energy by varying 
the electrode-electrode distance. We remark that in the 
optimized geometry, the plane of the molecule is approx- 
imately perpendicular to the plane of the two graphene 
sheets. This minimizes the coupling of ir states of the 
molecule with those of graphene, which can be advanta- 
geous in certain cases. In fact, it would be interesting to 
tune the orientation of the molecular plane with respect 
to that of graphene, by linking the molecule to more than 
one carbon of the graphene. In this way, one could ad- 
just the hybridization of the molecular levels by changing 
the relative orientation of the molecule with respect to 
the graphene plane. The system may thus be thought 
of comprising of three separate parts: the left and right 
semi-infinite graphene leads, and an extended-molecule 
region. For each electrode, six layers of carbon atoms 
are included in the extended molecule region in order 
to screen the perturbation from the central region. With 
this, the final geometry with the contacts is relaxed again. 




FIG. 1. Structures of the closed and open forms of the two 
diarylethene derivatives studied in this work. Here, R is the 
functional group that links the molecule to graphene, which is 
-CH 2 - in the first case (Al and Bl) and -CO-NH- (amide 
linking) in the second case (A2 and B2) 



The graphene electrodes are constructed so that the pe- 
riodic replicas of the contact region are separated by four 
unit cells (i.e. 7.436 A), and the interaction among them 
is negligible. Our electronic structure calculations are 
carried out using the DFT self-consistent pseudopoten- 
tial method implemented in SIESTA package [3(J. The 
exchange-correlation effect and electron-ion interaction 
is described by the Perdew-Burke-Ernzcrhof generalized 
gradient approximation and the norm-conserving pseu- 
dopotcntial in the fully nonlocal form, respectively. A 
double-C plus polarization basis set is used to describe 
the localized atomic orbitals and an energy cutoff for 
real space mesh size is set to be 200 Ry. The ionic 
positions of carbon atoms are the ideal positions of a 
previously relaxed graphene sheet, while the hydrogen 
positions are relaxed for each junction. We performed 
all relaxations with a force tolerance of 0.04 eV/A. For 
the self-consistent calculations, 12 K-points are used in 
the transversal direction (that of the electrode's edge), 
while for the calculations of the transmission function 
60 K-points are taken into account. The vacuum por- 
tions between two neighboring layers are set to be 15 A. 
The transmission function of the system is calculated 
by the TranSiesta code [3l| which combines the NEGF 
technique with DFT. The results are similar for the four 
molecules shown in Fig. [T] and we choose to discuss those 
of Bl in depth. 

In the upper panels of Fig. [2] we present the total DOS 
(solid line) of the extended system and the projected den- 
sity of states (PDOS) (shaded area) onto the molecular 
region, as function of the electron energy, the left and the 
right one for the closed and open isomer, respectively. 
The transmission function dependence on the injected 
electron energy is shown in logarithmic scale in the lower 
panels of Fig. [2] The total DOS is reminescent of that 
of the electrodes, since near Ep it displays a linear dis- 
persion following that of an ideal, infinite graphene sheet. 
Of course, very close to Ep the dispersion is not perfectly 
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FIG. 2. PDOS and transmission function for the closed (left) 
and open (right) isomers of Bl diarylethene. Upper panel: 
DOS (line) and PDOS (shaded area) on the molecular region. 
Lower panels: zero-bias transmission functions. 



Closed - LUMO 




Open - LUMO 




FIG. 3. Real part of some representative wavefunctions of the 
junction with molecule Bl: the LUMO of the closed isomer 
(upper panel) , the LUMO of the open isomer (middle panel) , 
and one edge state for the open isomer (lower panel). 



linear due to the presence of a vacuum region between 
two graphene sheets. Furthemore we remark the pres- 
ence of two low peaks very near Ep, one slightly above 
and the other just below it. They are states mainly lo- 
calized at the molecule-electrode interface whose origin 
stems from the zig-zag termination edge states of the 
two graphene sheets interacting with the orbitals of the 
molecule in the coupled system. We point out that, as 
a consequence of the perturbation due to the junction, 
such states lie within the continuum graphene spectrum 
and they become resonant conducting states, which do 
not exist in SWCNT |27|. Apart from the edge states, 
all peaks in the PDOS refer to molecular orbitals broad- 



ened by the interaction with the electrodes and a cor- 
responding feature at the same energy appears in the 
transmission function. We concentrate on the two main 
peaks of conductance for each isomer that appear on the 
left and on the right side of Ep. They are due to the 
hybridization of the HOMO and LUMO states of the 
molecule with the graphene substrate. For the closed 
isomer such states are closer to Ep, and their intensity 
is lower, if compared with that of other molecular states, 
since the DOS of the underlying graphene sheets tends to 
zero as the energy approaches Ep. However, considering 
small enough energy deviations from Ep, the transmis- 
sion function is always larger for the closed isomers than 
for the open ones. The lower conductance of the open 
isomers is due to their distorted geometry: the breaking 
of the central C— C bond due to the closed-open transi- 
tion affects the 7r-conjugation of the structure reducing it 
[3^ |. Hence, the electrons are less delocalized along the 
molecular backbone reducing transmission. In Fig. [3] we 
display the real part of the wavefunction for three rep- 
resentative states of the system, namely the LUMO of 
the closed isomer, the LUMO of the open isomer, and 
one edge state of the open molecule. We observe that 
the former LUMO is more widely spread than the lat- 
ter one, reflecting its more extended ir conjugation. The 
plotted edge state shows a rapid decay inside the molec- 
ular region. The other edge state just mirrors it on the 
opposite side of the junction. Comparing the transmis- 
sion functions of diarylethene with different connecting 
groups (Bl and B2), we noticed that the presence of the 
amide end-groups in the latter molecule produces a slight 
shift of the peaks toward higher energies. This results in 
a better alignment of the HOMO of the closed isomer 
with Ep so, once a voltage is applied, a larger current 
may be obtained in this case. But, as the HOMO state 
is shifted nearer to Ep, its conductance intensity is low- 
ered because the DOS's of the electrodes diminish. Apart 
from this shift, the structure of the transmission function 
is not appreciably affected by the change of connecting 
groups. 

We now discuss the changes in the conductance for 
molecule Bl when the system is driven out of equilibrium 
by the application of a bias V. In Fig. [4] the transmission 
functions are reported for the closed (left panel) and open 
isomer (right panel) by increasing the voltage from the 
bottom pictures. For the closed molecule a bias of 1 V 
allows one to include the contribution of the HOMO in 
the admissible energy window (vertical lines), while a 
bias of 2 V is needed to contain the LUMO. We also 
observe two features from the edge resonances near Ep. 
In fact, as a finite bias is applied, edge states play a 
role, though small, in the transmission curve. A non 
negligible conductance asks for a bias V = 4 V for the 
open isomer. In both cases the HOMO-LUMO energy 
difference is unaffected by the applied voltage. 

In solution, diarylethene based molecules can be 
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FIG. 4. Transmission function of the closed (left) and open 
(right) isomers of the Bl diarylethene derivative at different 
applied biases (from bottom to top: IV, 2 V, 4V). The bias 
window is delimited by two vertical lines. 



switched between the two states by photochemical 
means. When used as a bridge between two Au elec- 
trodes, the switching occurs, usually in one direction 
only, due to quenching by the electrodes. For CNT elec- 
trodes too, the switching from open to closed isomers 
occurs through light, but the reverse does not [lq ]. We 
suggest the use of graphene as the electrode as (1) it has 
a lower density of states near the Fermi level in com- 
parison with Au leading to lesser probability for quench- 
ing and (2) the orientation of the molecule with respect 
to the graphene sheet can be changed chemically, tun- 
ing the interaction between the two components. In 
the polymeric state, the change from closed state to the 
open state can be induced not only by light, but also 
by electron/hole injection obtained by the application 
of a potential difference Q. Consequently one may an- 
ticipate that in the case of diarylethene molecules con- 
nected to graphene, the switching from open to closed 
may be induced by light and the reverse by the appli- 
cation of a large potential difference. Looking at the 
positions of the LUMO/HOMO orbitals energies in the 
transmission curves given in Fig. [4j it seems possible to 
populate/depopulate them to a sizeable extent by the ap- 
plication of about 2 and 4 volts for the closed and open 
configurations, respectively. Thus one has the possibility 
of inducing switching in both directions, by the applica- 
tion of such a potential difference for a short period of 
time, and hence light may be avoided altogether. 

In conclusion we have calculated the conductance of 
the two isomers of diarylethene-based molecules and an- 
alyzed the role of the graphene leads. We have pointed 
out the advantages of using graphene, namely a re- 
duced quenching of the photoexcited state compared 
to gold electrodes and the possibility to tune the reso- 
nant HOMO and LUMO hybridization by orientating the 



molecule with respect to the leads plane. We also stress 
that, in the case of graphene, the switch can be operated 
electrically by the application of short pulses, making it 
potentially very useful, in comparison with other elec- 
trodes. 
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